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ABSTRACT 

Non-linear force-free fields are the most general case of force-free fields, but 
the hardest to model as well. There are numerous methods of computing such 
fields by extrapolating vector magnetograms from the photosphere, but very few 
attempts have so far made quantitative use of coronal morphology. We present 
a method to make such quantitative use of X-Ray and EUV images of coronal 
loops. Each individual loop is fit to a field line of a linear force-free field, allowing 
the estimation of the field line's twist, three-dimensional geometry and the field 
strength along it. 

We assess the validity of such a reconstruction since the actual corona is 
probably not a linear force-free field and that the superposition of linear force- 
free fields is generally not itself a force-free field. To do so, we perform a series of 



tests on non-linear force- free fields, described in iLow &: Loul (119901 ). For model 
loops we project field lines onto the photosphere. We compare several results 
of the method with the original field, in particular the three-dimensional loop 
shapes, local twist (coronal a), distribution of twist in the model photosphere 
and strength of the magnetic field. We find that, (i) for these trial fields, the 
method reconstructs twist with mean absolute deviation of at most 15% of the 
range of photospheric twist, (ii) that heights of the loops are reconstructed with 
mean absolute deviation of at most 5% of the range of trial heights and (iii) that 
the magnitude of non-potential contribution to photospheric field is reconstructed 
with mean absolute deviation of at most 10% of the maximal value. 



1. Introduction 

Most active region coronal magnetic fields are believed to be in a force-free state, 

V X B = a(r)B, 



(1) 
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where a is a scalar of proportionality (e.g., iNakagawa et al.lll97ll ). It turns out that a is 
closely related to the local twist of magnetic field lines. For example, in a cylindrical 



uniformly-twisted flux tube, field lines twist about the axis by an angle Q 
axial distance L (Aschwanden, 2006). 



\aL over 



If a varies in space, a solution of Eq. ([T]) is called a non-linear force-free field, as it 
solves a non-linear sy stem of equations for different components of B and the scalar a. 
Demoulin et al.l (119971 ) described basic problems arising when trying to solve these equations. 
In particular, the existence and uniqueness of a solution is not entirely clear. Another 
difficulty is the fact that the only source of boundary conditions available at the moment are 
vector magnetograms measuring the field within the non-force-free photospheric layer. 

A particular case of a force-free field, called a linear force-free field, or a constant-a field, 
occurs when Va = 0. In this case, using V ■ B = 0, Eq. ([T]) is transformed to a Helmholtz 
equation for B. This is much easier to solve and the conditions for existence and uniqueness 
of solution are known. Nor does the solution require vector magnetogram data, but only a 
line-of-sight magnetogram and a value of the constant a. This boundary condition is affected 
less by the fact that the photosphere is pr obably not force-free. There are many methods 



of solving for linear force- f r ee fields (e.g., iNakagawa &: Raadu 



Altschuler fc Newkirklll969l : iLothian fc Browning 



19951: 



1972 



Chiu fc Hilton 1977 



Alissandrakia Il981[). In par t icular 



in this paper we will use the Green's function method, described in IChiu &: HiltonI (119771 ). 
as it does not place any restrictions on a and it generates a field over an entire half-space, 
without boundaries. 

While they are simpler to generate, linear force-free fields have proven insufficient to 
model complex geometries of the solar corona. Observations of curvature of Ha structures, 
visual studies of twist in coronal loops, and estimations of local twist at the photospheric 
level via vector magnetograms rev eal active regions wi t h spatially varying twist, and even 
varying sign of twist (for example, iBurnette et al.l 12004 INakagawa fc Raadulll972l ). In light 
of this any constant-a approximation would appear to be, strictly speaking, incorrect. Over 
the past decade there have been many attempts to perform extrapolations of non-linear force- 
free magnetic field into the corona and to assess th e quality of the extrapolation by comparing 



Schrijver et al.ll2008l : iDeRosa et al 



2009 



lines of the resulting field to coronal loops (e.g. 
and references therein). 

In the current paper we apply a completely different approach. We use the visible shapes 
of coronal loops to infer the twist of the magnetic field. Instead of measuring twist in the 
photosphere, where Eq. ([T]) is not appropriate, we perform measurements in the region of 
interest, in the low-/? force-free solar corona. The method thus relies solely on EUV or X-Ray 
images of coronal loops and on line-of-sight magnetograms. 
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The basic idea is to try to approximate every visible coronal loop with a field line from 
a linear force-free field, and allow a to be different for every loop. Even for non-linear force- 
free fields, a must be constant along each field lin^. If a changes smoothly then it would 
be reasonable to expect that a is nearly constant in the vicinity of a given field line. Of 
course, a superposition of constant-a fields would not in general be a force-free field and at 
first sight such a method could not be expected to yield meaningful results. In the next few 
paragraphs we argue that such a method might work under certain circumstances relevant 
to the solar corona. Within the core of this work we support the hypothesis with tests first 
using analytic non-linear force-free fields and then with solar data. 

Our method is similar to the ones proposed by Green et al (2002) and Lim et al (2007), 
however, with several important advantages. First, it does not require the full length of a 
coronal loop to be visible for a successful reconstruction. Second, it does not require either 
of the footpoints to be visible. Third, it allows the user to draw a smooth curve (a Bezier 
spline) interactively on top of the loop, rather than selecting a few points along the loop. 
This maximizes the amount of information taken from the coronal image. The fit itself is 
similar to the one used by van Ballegooijen (2004); but while van Ballegooijen (2004) fits 
loops with lines of a particular non-linear force-free field model, we fit loops with lines of 
many different linear force-free fields, choosing the best a. 

Consider an imaginary example of two dipoles far apart compared to their sizes. Suppose 
that they constitute a non-linear equilibrium, having different twist, possibly of opposite 
signs. Suppose, however, that within each dipole the twist is more or less constant. In 
such a scenario there would probably be some transition region between the dipoles where 
a changed significantly. Provided the dipoles are far apart we may claim that in the vicinity 
of the footpoints of one of them the current of the other would not significantly perturb the 
field, and in the close vicinity of each of them the field would be nearly a constant a field. 

To support this reasoning we note that the dipolar term of a magnetic field drops as 
and thus the effect of a distant dipole is in general not very large compared to the nearby 
dipole. Indeed, this is why studying a magnetic field of an isolated region of the corona is 
at all meaningful. 

We therefore argue that a non-linear force-free field could be considered to be linear 
in the regions of slowly changing a (in some sense of the term). Thus the geometry of a 
field in isolated regions of slowly- changing a might be approximated by the geometry of a 
constant-a field. 



^This result is obtained by taking the divergence of both sides of Eq. ([T]) and using V • B = and the 
identity V • (V x B) = 0, so that B • Va = 0. 
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What is the limit of applicability of such an assumption? It is quite clear that it could 
work well for an isolated uniformly-twisted active region. Is it possible to pick a field line 
in an active region, and to suppose that the field's geometry is not significantly different 
from that of a linear field in a close vicinity of this field line? We herein conduct several 
experiments on both synthetic and real data which provide evidence that at least in certain 
cases of interest such an assumption is reasonable. 

The paper is organized in the following manner. In Section [2] we explicitly define the 
function to be minimized in order to obtain a best fit, a "distance" between two curves, d. 
In Section [3] we describe the minimization process, varying a and the line-of-sight coordinate 
/i, and report that it indeed works in the obvious case, where the loop is a field line from 
an actual linear force-free field. In Section H] we describe the typical features of the function 
d{a,h) and attempt to explain their appearance. In Section Owe present the r esults of ap- 
plying this procedure to several analytic non-linear force-free fields, described by lLow fc Lou 
(Il990l ). We also present an additional step proven necessary for the best fit procedure. This 
step amounts to minimizing d{a, h) in a very specific region of (a, h) parameter space. We 
demonstrate that this step significantly improves the results for strongly twisted fields. In 
Section E] we demonstrate the same method applied to real data: line-of-sight magnetograms 
from SOHO/MDI and coronal images from Hinode/XRT. In Section [7] we discuss the results 
and their possible use in studying coronal magnetic fields. 
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2. The distance between two lines 

In order to compare two curves, we seek a function quantifying tlie discrepancy between 
the curves. The ideal function would yield similar results to those obtained by visual com- 
parison. It would be large when a human observer would consider the two lines to be far 
apart or unlike one another and small when a human observer would consider them to be 
similar and close to one another. 

We use for this purpose a function first introduced by Green et al (2002) and later used 
by Lim et al (2007). We apply it, however, to a different set of objects. While Green et 
al (2002) and Lim et al (2007) compare a few points along the loop to a set of field lines 
traced from the photosphere, at the presumed location of the loop 's footpoints, we compare a 
smooth curve, chosen to visually match the loop, to a set of field lines traced from different 
loc ations along the l i ne-of- siqht at some point along the curve. The same method was used 



by Ivan Ballegooijenl (j2004j ). Therefore, unlike method of Lim et al (2007), our method does 
not require knowledge of the footpoints. In fact, it will work with even a small portion of a 
loop. 

The discrepancy function is defined between two smooth curves in a plane, Li = 
{xi(/), ?/i(/), < / < Li} and L2 = {x2{l),y2i})i^ < I < L2}. For every point / on Li it 
is possible to define a minimal distance between that point Li(/) and L2 in the classical 
sense: the smallest of the distances between the point Li(/) and every point along L2. It 
could also be defined as the length of the shortest perpendicular from the point Li (/) to the 
curve L2, given sufficient smoothness. We will refer to this distance as (5(Li(Z), L2). The 
discrepancy between the two curves is the average of 5{l) over curve Li, 

Li 

d{UM) = ^J 5(Li(/),L2)c^L (2) 


A numerical scheme to compute this integral is quite simple. Sample line Li into ui 
segments with equal lengths A/i = Li/ui. Provided the segments are small compared to the 
local radius of curvature, for every point on Li, 

5{U{k),U) = min (^^{xi{U)-X2Y + {yi{k)-y2Y\ {X2,y2) G L2) . (3) 
The discrepancy integral can then be approximated by the sum 

A/ 

d{UM) = -r^Y.^{Mk)M), (4) 

J^i — 



i=0 
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representing the mean distance between points of one curve and the whole second curve; it 
has the units of length. 

The trivial properties of the discrepancy function (i(Li,L2) are, first, that it is non- 
negative and second, it is non-commmutative, meaning d{hi, L2) is different in general from 
(i(L2,Li), as illustrated in Fig. [H 



d{L., LJ = 0.447 d{U, L,) = C.705 




0246 0246 



Fig. 1. — The discrepancy between the two lines. The left panel illustrates the calculation of 
d(Li, L2) and the right shows (i(L2, Li). For each point on the first curve (the first argument of d), 
denoted by a triangle, one finds the closest point on the second curve (the second argument). The 
closest distance is shown as a dashed line. The net discrepancy is the average of all such distances. 
Comparison of the two panels illustrates that the discrepancy is non-commutative. 



3. a-h-fit 

For a visible coronal loop it is possible to construct a smooth two-dimensional curve Lq 
in the plane of the sky visually approximating the loop, or some portion of the loop. The 
loop is really a three-dimensional structure and for every point on Lq, the third coordinate, 
i.e. along the line of sight (LOS), is unknown. If the loop is at disk center, then the LOS 
coordinate is the height above the photosphere. For simplicity in notation, we will thereafter 
refer to this coordinate as "height" , denoted h, even when it is not vertical. 

The main idea of what we call an h-fit is to choose a point Iq on the loop Lq and 
prescribe a certain height. Then, if the magnetic model is known, trace a field line from the 
three-dimensional location {Lqx{Io) , L^yil^) , h) and compare its plane of the sky projection 
L(/i) to the original loop by calculating d{h) = (i(Lo, L(/i)). Finally, we vary h to find the 
minimum of d{h). 
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To illustrate this method we construct a synthetic magnetogram and generate a magnetic 
field (the potential field of a magnetic quadrupole in half-space Z+). As a model of the 'loop' 
we take an actual field line, projected onto the x-y plane. We then take the mid-point of the 
projected loop, Lo(^Lo), and trace field lines at different heights (see Fig. [2]). Fig. [3] shows 
the function d{h) with one minimum at the actual height, to within one step of the h search. 

Suppose now that the coronal magnetic field is not known, but belongs to a known 
family of magnetic fields, described by certain parameters. It is then possible to do fitting 
not only in height, but also in the space of these magnetic field parameters. For example, if 
there is a reason to believe that the actual field is constant-a , but with unknown a, then 
the discrepancy is a function of both a and h, d = d{a, h) and the minimization must be 
done over (a, h) space. 

To illustrate such a minimization we perform the following experiment. From the 
quadrupolar magnetogram of Fig. [2] we extrapolated a constant-a field into the corona (2'+). 
We produced several such fields with different values of a. For each of these fields we selected 
several magnetic field lines at random, which we projected onto x-y plane. We then treated 
these projections as synthetic loops and performed an a — h fit. Provided the method works, 
the best-match a should have one-to-one correspondence with the real a for each field line. 



The constant-a fields were generated using Green's function (jChiu &: HiltonI 119771 ) for 



a field in a half-space. This has the advantage that it places no limitatio ns on a, whereas 



Fourier methods using periodic images separated by L, require |a| < ^ ( iNakagawa et al. 



19711 ). In Section [5] we demonstrate that our method can reconstruct values of |a| up to ^, 
where the height h of a point (not necessarily the highest one) along the loop could be much 
less than the linear size of the computational domain L. For the method we propose, to 
perform a — h fit within the full region of interest, including h < L/2, a values larger than 
the maximal allowed by the Fourier method are needed. 

The results of the fit show that there is indeed a strong correlation between best-match 
a and the real a. However, in some cases (46 points out of 689, about 7% of all field lines) 
the fit seems to be off by more than one step of a. We attribute these errors of the fit to 
several factors. One factor is the gridded search algorithm whereby we calculate d{a, h) for 
each point on a grid with fixed steps in both a and h (it is clear that a better algorithm could 
be implemented, however, in this paper we concentrate on the theoretical possibility of the 
method, rather than on programming tasks). A second is that the fit is poorly constrained 
when a field line's shape is hidden by the projection. Finally, there are numerical errors 
associated with numerical integration of a field line from a field represented only on discrete 
grid points. 
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4. Shape of d{a, h) in the parameter space 

Fig. m sliows the function d{a, h) for one of the loops from the experiment described 
above. This function has vaUeys (dark) in the shapes of hyperbolae, located at or between 
the hyperbolae h = nir/a, n = ±1, ±2, ... After examining parameter spaces of many field 
lines we have concluded that there is one and only one "branch" of local minima in each 
nix /a < h < {n + l)TT/a, except for n = ±1; there are usually two or more "branches" in 
—7i/a < h < n/a. 

The foregoing behavior can be explained by the Green's function in the far field. Far from 
the photospheric fiux concentration at fq the field is dominated by the monopole moment, 
B(r) (X G(r — rp). For a field restricted to half spa c e, Gr een's function G is given in 
Chiu &: HiltonI (119771 ) and later in iLothian fc Browningl (Il995l ). It depends on sm{ujh + 0o) 



and cos{ujh + (po), where u oc a and (f)o oc a^/{x — xq^ + {y — UoY- If one changes both 
a and h in such a manner that ah = const, that is, remaining on the same hyperbola in 
{a, h) parameter space, the sm{ujh + 0o) and cos{ujh + 0o) terms remain constant. 

This observation about the parameter space suggests a useful heuristic restriction to the 
search. Tracing field lines from {xo,yo, h) for increasing values of h (provided a 7^ 0) causes 
the angle between the field line and L at (xq, Vq, h) to increase or decrease monotonically. 
The cosine of this angle, 

B(a;o,yo, /i,a) ■ L(a;o, 2/o) 



cos(^) 



IBIILI 



will be a local maximum near the local minima of d{a, h), but not exactly at the same place. 
Evaluating the magnetic field at one point is, of course, much faster, than tracing a whole 
field line. We found that without loss of any information about local minima of d, we may 
restrict the search to only those {a, h), for which | cos(^)| 

It is clear why the local minima of d{a, h) are co-aligned with maxima in | cos(6')|. Since 
L is a smooth curve, there will be a vicinity of (xq, 1/0), where L is, to first order, a straight 
line. The same is true for field lines of a force-free field: the field lines are smooth curves, so 
the field line traced from {xo,yo, h) in a close vicinity of this point is also to the first order 
a straight line. Suppose |cos(6')| = 1, i.e., B(xo,i/o;^) is parallel to L(xo,?/o)- Then in the 
neighborhood of (a;o,yo) the two curves would be exactly the same, and d would be close to 
zero, if averaged only in that vicinity. Farther from (xo,yo) the two may differ significantly 
resulting in non-zero d over the whole length of the loop. 

If the line-of-sight angle is such that the loop's projection is nearly a straight line, then 
there will be many field lines, with different a, that are high enough and long enough to 
appear nearly straight along all the length of the projection. In this case the fit may give 
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poor results. The upshot is that even though one does not need the whole loop to perform 
the a-h fit, the projection of the visible portion should not be "too straight". We develop a 
more quantitative measure for this criterion below. 

5. a-h fit: applied to Low & Lou field 

We next test the a-h fit on a set of non-linear force-free fields from Low & Lou 
(1990). Each of these can be viewed as the field of a singular point source placed below 
the photosphere and inclined. The field is specified by parameters a (related to character- 
istic range of field's a), I (depth of the source under the photosphere, we used / = 0.3 for 
all fields), $ (orientation of the source, we used $ = 7r/2 for all fields) and n (for explicit 
derivation and definitions, see the Appendix). 

The experiments were conducted as follows. We generated several Low & Lou fields with 
different parameters. For each field, we traced a few hundred field lines, projected them on 
the x-y plane and used them as synthetic loops. Then for each such loop we conducted an 
a-h fit (by gridded search), with values of a being within the range of Low & Lou's field 
photospheric a. 

We found that for a "dipolar" field (in the sense of it having two distinct polarities, see 
Fig. [5]) with n = 1, a = 0.02, the values resulting from constant a fits do indeed correlate 
with the real values of the field lines, as shown in Fig. [6] (left plot). 

For this dipolar field the d{a, h) plots for most loops had one distinct valley of local 
minima (horizontal or nearly horizontal), and hints of other valleys at larger a or h. Another 
notable feature of the low valley is that it tends to cross the a = line rather than approach 
it asymptotically. The parameter space for one of these loops is shown in Fig. [7] (left plot). 
Notable in that plot is the global minimum was always in the lowest nearly-horizontal valley. 
The "true location" (from the original Low & Lou field line) was also within that valley; 
however, it is sometimes offset with respect to the global minimum. In general, a of the 
global minima are correlated with a of the original Low & Lou field lines, as shown in Fig. [6] 
(left plot). Finally, Fig. [7] (right plot) shows that the constant a field line of the global 
minimum seems to approximate the original Low & Lou field line quite well, although a 
tendency to under-estimate a is evident, and still clearer in the histogram of Fig. O (right 
plot). 

For notational convenience, we hereafter refer to lines of Low & Lou fields as real field 
lines and to their best-fits of constant a field as found field lines. We will also use (a^eaz, hreai) 
to denote the parameters of the real field line (recall: h is a. height in the midpoint of the 
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line's projection into the photosphere). Similarly, we will use [a found, h found) notation to 
refer to the parameters of the found field line. 

Best-fits are potentially useful in reconstructing the photospheric distribution of a. We 
constructed a photospheric map of a by assigning the coronal value to the footpoints of the 
reconstructed field lines. A full map requires a smoothing, averaging or interpolation, to 
assign a to photopsheric points around the footpoints of observed loops. To illustrate this 
possibility we did a robust reconstruction with bicubic spline interpolation (see, for example. 



Press et al.lll986l ). shown in Fig. [8l The fit that we did is simple and robust, nevertheless it 
is able to reconstruct the general shape of the actual distribution of a in the Low and Lou 
field. 

Another measure of the quality of the fit is its reconstruction of magnetic field B, for 
example, at z = 0. We utilize the form Bjj = Bpot+B„p, where Bjj is the full force-free field, 
Bpof is the potential field with the same normal component at the boundary, and B„p is a 
"current contribution" — a non-potential force-free field with B„p- n|av = 0. Note that Bpot 
is the same for ^l&cL and Brecon (reconstructed), since it is uniquely defined by the volume 
and by the Dirichlet boundary conditions. For "weakly non-potential" field |Bpo(| ^ |B„p|; 
this is true of some of our cases. Rather than comparing ^l&cL to Brecon, we compare their 
"current contribution" terms, normalized by the potential field: IB^^^ — Bpot|/|Bpot| to 
l^recon — Bpot\ / \Bpot\- The histogram for z G [0,2] pix is shown in Fig. [9] (left). To make 
it, we evaluated ^l&zL, Brecon and Bpot along the found field lines. It seems that for most of 
points the two fields were nearly identical, suggesting the accuracy of the reconstruction. 

From both Fig. [6]and Fig. [9]it seems that the reconstruction does a better job for smaller 
a and for weaker B„p, than for larger a and for stronger B„p. This and the reasoning from 
the previous section suggest that a-h fit might not work for strongly twisted, or maybe strongly 
non-linear, fields. We tried to determine the range of a for which the fit would yield reliable 
results. For that, we generated several more Low & Lou fields, this time quadrupolar (in 
the sense of it having three polarities, like the field of a point quadrupole): we kept n = 2 
and gradually increased a over the following values: a G [0.05,0.1885,0.3,0.6,1.0,1.5,2.0]; 
in addition we computed a field for = 3, a = 0.4. We generated both signed and unsigned 
Low & Lou fields; both have identical photospheric B^, but the first one has a > and 
a < 0, while the second one has only a > (see Appendix). This was done in order to relate 
the errors of the fit with the "non-linearity" , that is, with how much a changes over a fixed 
length. 

Typical parameter spaces for those fields are shown in Fig. (TU] and Fig. [TTJ We found 
that for weakly twisted fields (that is, a = 0.05, a = 0.1885 and marginally a = 0.3, for which 
|Brecon — Bpot|/|Bpot| is at most 0.15, 0.3 and 0.4 respectively within the computational box 
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close to the photosphere) a typical d{a, h) parameter space has one valley of local minima, 
with the same characteristics as the parameter space for n = 1.0, described earlier. It does 
not seem to approach a = or h = asymptotically like a hyperbola would, but rather it 
crosses a = axis. The global minimum always lies on this valley. The point {areai, Keai) 
also lies on this valley. It seems that the more horizontal this valley is, the more offset could 
be the "real" and "found" points along the valley, so the more different could be areai and 
a found, Keai and h found always appear to be very close. 

For more strongly twisted fields (a > 0.3) the parameter space within the range of a and 
h of the real field reveals more valleys, although the larger-scale behaviour seems to follow 
the analytic Green's function behavior shown in Fig. HI That is, d{a, h) seems to have valleys 
of local minima that look like hyperbolae and are located at or in between the hyperbolae 
ah = ±7r, ±27r, .... Except for \ah\ < tt, there is one and only one valley in between every 
two hyperbolae mr < ah < (n + 1)tt, n = ±1, ±2, .... Within \ah\ < tt there are usually two 
or more valleys, and one of them is usually "non-hyperbolic" in the sense described above. 
In these more strongly twisted cases we observed that the global minimum could be in one 
of the "higher" valleys. It seems that the field line corresponding to [a, h) of the global 
minimum is much longer than the "loop" (line of Low & Lou field) and morphologically is 
quite different. Its smaller d results from a small portion of the long line coinciding with the 
"loop" . This happens especially often for loops that are "too straight" in some sense. This 
is qualitatively described in Section H] and quantitatively described further in the text. 

After inspecting a great number of these parameter space plots we have noticed that 
{areai, Keai) still tcuds to Correlate with the location of the "non-hyperbolic" valley. To prove 
this point, we conducted the following experiment. First we excluded loops that were "too 
straight" . Second, we chose as the best fit for each of the loops a local minimum on the non- 
hyperbolic valley, rather than the global miminum. The results of this two-step procedure 
are shown in Fig. [T31 The explicit description of the procedure is below. 

As a definition of "too straight" we adopted the ratio of sides of a box circumscribing 
the loop. The box is aligned with the least-square line fit to the loop, its length being the 
length of the loop along this line and its width being twice the maximal deviation. Based on 
visual examination we chose the minimum width-to-length ratio to be 0.05 for the "loop" to 
be eligible for the analysis. 

As for selection of the "non-hyperbolic" valley, we developed and followed an algorithm 
based on the shape of the parameter space. We have found that for Low & Lou fields this 
algorithm yields good results. First of all, for a given parameter space plot we identified 
several one-dimensional local minima for each column {ai,h). Then we manually select 
some of those local minima that belong to only one of the valleys and find a local minimum 
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of d{a, h) within this valley. 

For the selection of the valley, we followed these steps: 

0. Consider only the valleys for which \ah\ < vr. 

1. Is there one "non- hyperbolic" valley in this region? If yes, select the local minima 
within it. If definitely no, proceed to the next step. If not sure, discard this loop from 
consideration. If there are several local minima along this valley, select the one that 
has the lowest h. Example in Fig. [T21 top left. 

la. If the "non-hyperbolic valley" merges with a "hyperbolic- like" loop, select the local 
minima in the "non- hyperbolic" part. If unclear, discard this loop from the consider- 
ation. Example in Fig. [121 top right. 

2. Does this "non- hyperbolic" valley seem to change directions, possibly crossing a = 
more than once? If yes, select local minima on the lowest (smallest h) section of it. If 
definitely no, proceed to the next step. If not sure, discard this loop from consideration. 
Example in Fig. [121 middle left. 

3. Are there two "non-hyperbolic" valleys on either side of a = 0, and neither of them 
crosses a = line? If yes, select local minima on the one that extends to a bigger 
range of h. If definitely no, proceed to the next step. If not sure, discard this loop 
from consideration. (We found that such parameter space plots often happens for a 
"too straight" loop, and threshold of 2a/L = 0.05 seems to eliminate the majority of 
them. For the latter ones, areai seems to be on the higher-extending valley.) Example 
in Fig. [121 middle right. 

4. If there is no such special valley, among the "hyperbolic" valleys in \ah\ < vr there is a 
"lowest-order" one, that is, the one that has smallest \a\ for h — *■ oo. Is there enough 
of this loop presented? (I.e., that did not fall below the threshold on \cos9\ < n/2, 
as described in Section |H or that did not fall below any other threshold that was 
used, such as difference in length being too big, or the length of the field line being 
significantly smaller than the length of the loop, or the amount of self-crossovers of a 
field line being two large - we use the second and third thresholds, but not the first 
one.) If yes, select the lowest in h local minima on this valley. If definitely no or not 
sure, discard this loop from consideration. Example in Fig. [T^l bottom left. 

5. Hard to classify cases: discard from consideration. Example in Fig. [121 bottom right. 
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As shown in Fig. [131 for signed field with a = 1.5, the global minima selection does 
not work very well. The above-mentioned algorithm of selection of only "non-hyperbolic" 
minima works much better; it significantly improves the correlation of areai and a found for 
large a (and big ranges of a). We also tested this algorithm for when the loops belong to 
linear force- free fields and verified that it yields the correct results at least within the range 
\a\L < 5, which is far beyond the range of all Low & Lou fields studied in this paper. 

The results for all Low & Lou fields are summarized in Table [TJ The individual results 
are shown in Figs. [H] - [T71 This includes scatter plots of areai versus a found, Keai versus 
h found, the comparison of B// and Brecon and photospheric distributions of areai and a found- 
The magnetic fields are compared in the same manner as described in Fig. (9) Bff and Brecon 
are evaluated at the photospheric level for each of the reconstructed field lines, and a two- 
dimensional histogram is computed. The photospheric distributions of a found are plotted 
in the same color table and with the same contours as areai and are obtained in the same 
manner as described in Fig. [8j afound{x,y, z = 0) is collected from all reconstructed field 
lines; the resulting set of points is used for two-dimensional spline interpolation. 

We draw several conclusions based on the results of this analysis. First, at least for 
some range of a, field lines of Low & Lou fields could indeed be approximated with the field 
lines of constant a fields of similar a and h. The reconstructed photospheric distribution of 
a seems to recover the general shape of the original field. Amazingly, it is also able to recover 
the area of the strongest gradients of a. Second, the height of the loops is reconstructed 
very well for the fields with a small range of a and less well for the fields with a greater 
range of a (see correlation coefficients and errors in Table [H note that for unsigned fields 
the range of a is about half the range in signed fields). Third, this method is also capable 
of reconstructing the magnetic field, at least near the photosphere. 

We summarize all the results on two plots in Fig. [TSl For each Low & Lou field, we 
looped through a found and measured the mean and standard deviation of areai, and plotted 
{areai) ± o" vcrsus a found- Wc did the same for h. It seems that the method systematically 
underestimates a by a small amount and it sometimes overestimates /i by a small amount. 
The least-squares line fit of the mean values, including standard deviation, gives an estimate 

(areai) OC 1.23afound and (hreal) OC 0.79hfound- 
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Side view 
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Fig. 2. — The synthetic example used to illustrate the h-fitting routine. Top: The synthetic 
magnetogram (in gray-scale) was used to generate the potential field in One field line was 
selected to represent the synthetic loop (in white; its starting point is shown as the triangle). It 
was projected onto x-y plane and compared with field lines, traced from the points with the same 
{xo,yo) (in this case {xo,yo) ~ (38,51)), but at different heights. These trial field lines are shown 
as dashed lines; they are traced for every pixel of the column {xo,yo,z e [l,nz — !])■ To make 
the plot clearer, only every second one is drawn. The point (xo,yo) was chosen as the mid-point 
of the projection of the "loop" (shown in diamond). Bottom: The same example, viewed in x-z 
projection. The initial "loop" is shown as gray, the traced field lines are shown as dashed lines, 
and their starting points are shown as diamonds. The thick dashed line shows the best h-fit. 
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Fig. 3. — The function d{h) has its minumum at a height h ^ 17 within the numerical error of 
the real height of the field line it is modeling. Note that for /i < 11 pix and /i > 93 pix the length 
of the projection of the traced lines is smaller than the length of the "coronal loop". Such lines 
could be automatically discarded from consideration, since the "loop" should be a part of a field 
line, and the length of a part of a curve cannot be greater than the length of the whole curve. We 
discard them by making d{h) artificially large if the length of Lq is smaller than the length of Li. 
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d(a, h)<0-2 for quodrupole field, loop #300 





Fig. 4. — (left) - The parameter space of d{a,h). The yellow asterisk shows the location of the 
"real" {a,h) of the field line and the cyan triangle shows the "best-fit" {a,h). The blue dotted 
lines are hyperbolae h = nvr/a, n = 0,ibl,ib2, ... For this field line, as for nearly all of them, the 
"real" minimum lies in the "branch" of local minima, that is within —n/a < h < n/a. We sped 
up the computation significantly by computing d only for (a,/i), for which at the initial point 
(xo,?/o) the magnetic field makes a relatively small angle with the normal of the loop: | cos(0)| = 
B • / /|B||Z| < V2/2. Our study shows that for most of the "loops" the local minima lie within this 
range of | cos(0)|. (right) - The parameter space of the Green's function for constant a field in half- 
space. The function being plotted is | cos(^)| = Gy{r — tq, a) / ^ G^{r — tq, a) + G^{r — tq, a) , 

where tq = (0,0,0) and ri = (10 cos(320°), 10sin(320°), h). The red contours are | cos(6')| = \/2/2, 
the blue dashes contours are h = nvr/a, n = 0, ±1, ±2, ... 
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Fig. 5. — (left) - A synthetic magnetogram of a non-linear nearly-dipolar field at the photosphere, 
with synthetic "loops" - field lines - projected onto x-y plane. For each of those lines we tried to 
approximate it with a line of a constant a field, recording a found and comparing it with the real 
a of the original field line, (right) - Best-fit field lines of constant a fields, each line belonging to 
a different constant a field. Hereafter all lengths are given in the units of Lq - a characteristic 
separation distance between two polarities, we calculated it as the distance between the pixels with 
maximal and minimal magnetic fields. In this case, Lq = 6 pix. 



Low & Lou field, n=1, a = 0.0200 
a^g^i versus df^^^^ scatter plot 
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Fig. 6. — Results for the data, shown in Fig. [5l On the left is a scatter plot of aj-eai vs. a found- 
The correlation is evident. On the right is the histogram of the error areai ~ found- 
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Fig. 7. — (left) - A typical d{a,h) parameter space for one of tlie field lines shown in Fig. O 
There is one nearly horizontal valley of local minima and a clue to other possible valleys for larger 
{a,h). White diamonds show the local minima in individual columns. Cyan triangle shows the 
location of the global minimum. A yellow triangle shows "real" (q, h) of the Low & Lou field line. 
Blue dashed lines (barely visible on this plot, but more evident on the other plots of this kind) are 
hyperbolae ah = nvr, n = 0, ±1,±2, ... (right) - Best-fit for the same field line (cyan), field of the 
constant a field, traced from the "real" values (yellow), and the the Low &: Lou field line (red). 
The difference is barely visible, however, the cyan line seems to match red line better than the 
yellow line. 
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REAL a 




FOUND a 






Fig. 8. — (left) - Photospheric distribution of a for the field from Fig. [5l The yellow contours 
are contours of a. (right) - the result of the reconstruction using a-h fit. The grayscale and 
contours on this reconstruction are identical to those on the real distribution. The red contours 
show the location of the footpoints of the found lines, that is, there is no information outside of 
these contours, and whatever is outside is shown solely for easy viewing. The only meaningful part 
is inside of the red contours. The result was extrapolated using thin plate splines fit into the set of 
footpoints with found a. This robust fit is sensitive to individual noisy points, and it is intended 
only to illustrate of the potential possibility of such reconstruction. Yet, with all these remarks, 
such robust fit is capable of reconstructing the principal shape of the distribution. 
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Fig. 9. — (left) - the comparison of reconstructed magnetic field, 'Brecon to the Low Sz Lou field, 
Bl&L) at the photospheric level {0 < h < 2 pix). We follow the "found" (best-fit) field lines 
and evaluate B;y^y (with different a for each field line), B^^^ and Bpot along them. By Brecon we 
mean the composition of B^jjj's for all loops (strictly speaking, the way it's constructed makes it in 
general not force- free and probably not even source-free, so it could hardly be called a magnetic field 
at all, rather, an approximation of the reconstructed field, evaluated along diff'erent field lines). As 
discussed in the text, most of the field is potential, so we compare the "non-potential" contributions 
only, normalized by the potential field. Here Bpot is a potential magnetic field, that is restricted 
to the upper half-space and has the same Dirichlet boundary conditions as Bi&i and the same as 
Brecon- It seems that the non-potential part of the magnetic field is reconstructed well, (right) - 
scatter plot of hreal vs. hfound- It seems, that h is found with much better confidence than a. 
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Low&Lou field, n = 2, a = 0.0500 I Loop § 




Low&Lou field, n = 2, = 0.0500 I Loop § 2 I 8.0 pix 




alpha [1/Lq] 

Low&Lou field, n = 2, = 0-3000 I Loop § 25 




if. 



-0.4 -0 2 0.0 0.2 

alpha [1/L^] 



Low&Lou field, n = 2, a = 0.6000 I Loop § 52 
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Low&Lou field, n = 2, □ = 0,1885 I Loop § 10 I L^= 8.Q pix 



Low&Lou fiold, n = 2, = 0,3000 I Loop § 25 I L„= 8,0 pix 




Low&Lou field, n = 2. = 0. 6000 I Loop fl 52 I L^= 7,0 pix 



Fig. 10. — (left column) - typical parameter spaces for Low & Lou fields with gradually increasing 
a. Here dashed blue lines are hyperbolae ah = 0, zLvr, ±27r, white diamonds show local minima 
in each column a = const, yellow triangle shows the location of {oij-eauWeai) and cyan triangle 
shows the location of the global minimum, that is, {oifoundih found)- (right column) - the original 
"loop" of each of those parameter spaces (red), the "global minimum" field line (cyan) that has 
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Qugdrupole field I Loop I 18 Qjodrupok tick! I Loop #18 1 L,-= 7.0 pix 
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Fig. 11. — The continuation of Fig. [10] for larger a, the notation is the same. Note that {oireaU Weal) 
is stiU at or near the "non-hyperbohc" vahey, while the global minimum could be at one of the 
other valleys. The field line, corresponding to the global minimum, is much longer than the loop 
and part of this field line happened to match the loop. 
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Fig. 12. — The illustrations of the algorithm that helps to select local minimum on "non- 
hyperbolic" valley. The notation is the same as in Fig. [TOl For description, please refer to the 
different options in the algorithm. 
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Low & Leu field, n = 2, a=1,5000 Low & Lou field, n = 2, a=1,5000 Loa' :V Loj -'icid, r =2, ;j=',5000 

areoi versus C(,(,^nd sca:ter plot a^eoi versus C(,(,^nd sca:ter plot a,^^! versus ^fa^na sootter plot 




Fig. 13. — (left) - a scatter plot of areal vs. ocfound for "signed" o = 1.5, n = 2 Low & Lou field. 
(middle) - same, but with "too straight" loops removed, (right) — same, but with "too straight" 
loops removed and with the minima selected only along "non-hyperbolic" valleys rather than the 
global minimum. 



- 25 - 




Fig. 14. — The results of the fit for "signed" Low & Lou fields: scatter plots of Oreai vs. a found 
(first column), hreal vs. h found (second column), histogram of |Bjj — Bpot|/|Bpot| (third column), 
real and reconstructed photospheric distributions of a (fourth and fifth columns respectively, the 
notation is same as in Fig. [8]). Note that as the range of a increases, the reconstructed photospheric 
distribution of a gets better, but the correlation of hreal and h found gets worse. 




Fig. 15. — The results of the fit for "signed" Low & Lou fields. The notation is the same as in 
Fig. [TH Note that the last row corresponds to an n = 3 field. 
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Fig. 16. — The results of the fit for "unsigned" Low & Lou fields. The notation is the same as in 

Fig.m 




Fig. 17. — The results of the fit for "unsigned" Low & Lou fields. The notation is the same as in 
Fig. [TH Note that the last row corresponds to an n = 3 field. 
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Table 1: A brief summary of the results, shown in Figs. [TUfTTl LAD fit and Spearman's rank order correlation 
were computed for areai{c( found)- LAD mean absolute deviation for a is given in percent of the range of a along the 
photosphere, for h in percent of the maximal height, and for |B — Bpot|/|Bpot| in percent of the maximal value along 
all measured field lines in the same height range that is used for the histograms. All statistics for the latter one is 
measured for |B — Bpoi|/|Bpof| > 0.05. The confidence of Spearman's rank-order correlation is bigger than 99.95% of 
a for signed n = 2, a = 0.05 and bigger than 99.999% for all other entries. 



-30- 



Compound data from all Low & Lou fields Compound data from all Low &c Lou fields 



J 0.0 
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Fig. 18. — (left) - Compound plot based on all Low & Lou fields measurements. For each field we 
went through all a found stnd measured mean and standard deviation of ctreal ■ Each such measure- 
ment corresponds to a point and error bar on this plot, (right) - same for h. The units for a are 
1/Lo and for h are Lq, where Lq is a separation distance between two polarities. 
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6. a-h-Fit Applied To Solar Data 

As a further test of the applicability of the method, we applied the algorithm described 
above to real solar data. We chose the unnumbered active region close to disk center on 2007 
February 10. This region, denoted SOL2007-02-10T11:22:L290 C101 according to t he Solar 
Object Locatoi§|, was observed by t he instruments on Hinode (Kosugi et al. 2007), partic- 
ularly the X-Ray Telescope (XRT, JColub et al.ll2007l : iKano et allboosh l As reported by 
McKenzie fc Canfieldl (120081 ). the region had a sigmoidal shape indicative of highly sheared 
coronal loops. Such regions are believed to comprise strongly twisted non-potential magnetic 
fields, and thus store significant amounts of energy, and commonly erupt in flares and/or 
coronal mass ejections. For this reason, coronal sigmoids present a very reasonable target 
structure for application of the proposed method of magnetic field modeling. The 2007 Febru- 
ary sigmoid is particularly useful because the high- resolution observations of Hinode/XRT 
(1 arcsecond per pixel) and the sigmoid's location near disk center facilitate tracing of the 
individual loops forming the sigmoid. We utilized the XRT image from 2007 February 10, 
11:22:06UT, to identify and trace 12 coronal loops. The magne tic models were gene rated 
from LOS magnetograms from the MDI instrument on SOHO (jScherrer et al.lll995l ). and 
then utilized for fittings to obtain the 3D shape and a of the coronal loops. 

To get the two-dimensional shape from XRT image we manually drew smooth curves 
(3-point cubic spline) over each of the loops. We then visually co-aligned XRT and MDI 
images and obtained (x, y) of the loops in MDI coordinates. No de-rotation was needed since 
MDI and XRT data were within one minute of each other. Since the region of interest was 
small enough, and close to disk center, we worked in a tangent plane approximation where 
the photospheric plane, 2; = 0, was taken to be the plane of the sky. 

To perform the fitting we extracted a region measuring 506" x 506" from the full-disk 
MDI magnetogram (257 x 257 pixels). In order to save computation time, we downsized the 
magnetogram by a factor of two in each dimensi on. Then we generate d constant a fields 
confined to half space (using Green's function from lChiu fc HiltonI (119771 )) in a 506" x 506" x 
200" computational box. We generated 41 different fields with —0.04 < a < 0.04 arcsec"^. 

The results are summarized in Table [2] and Fig. [191 Visually it seems that the fit did 
a good job for all but three loops. The parameter space plots, such as the one shown in 
Fig. [201 all looked like those for Low & Lou fields. For all of them we applied the proposed 
algorithm of selecting local minima on "non-hyperbolic" valley. 



Disregarding the three "unsuccessful" fits, it seems that a was of the order of 0.010 



^http://www. iac.es/proyecto/iau_divii/IAU-DivII/documents/target_naming_convention. html 
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0.015 arcsec"^ in the outer region of the sigmoid, and higher in the middle, exceeding 
0.02 arcsec"^. Using the solar radius 983.13" and recalling a typical bias of areai oc l.23a found 
(see Fig. [TS]), we estimate areai to be 1.7 — 2.6 x 10"®m~^ in the outer regions and over 
3.5 X 10~®m~^ in the core. These va l ue fall in the range typical of active region fields, such 



as those reported in iBurnette et al.l (120041 ): ±4 x 10~*m~^. The magnetic field strength 
along all the non-potential field lines was always within 50% of the strength of the potential 
field, evaluated along the same path; it was within 25% in at least half the cases. 

Assuming a typical separation distance L ^ 50" (see Fig. [19]) the values of a found was of 
the order of 0.5 — 0.75L~^ outside of the sigmoid and exceeded 1.0L~^ inside of the sigmoid. 
This is within the range of reconstructed a in the trial Low & Lou fields (see Fig. [18]) . 



7. Discussion 

In this work we have proposed a semi-automated method that, from a given two- 
dimensional EUV or X-Ray image of coronal loops and a line-of-sight magnetogram, re- 
constructs local twist and a three-dimensional shape of each loop. It tacitly assumes that 
coronal loops follow magnetic field lines. Our fitting matches the plane of sky projection of 
each loop to the projection of lines of linear force-free fields, traced from different heights 
along the line-of-sight and that have different twist in them. The method thus constructs a 
linear force-free field and one field line in it, that is the best match of an individual loop. 
The method is similar to the ones proposed by Green et al (2002) and Lim et al (2007), 
however, it holds several important advantages. First, it does not require a full length of the 
loop to be visible for a successful reconstruction. Secondly, it does not require any of the 
footpoints to be visible. Thirdly, it allows the user to draw a smooth curve (e.g., a Bezier 
spline) interactively on top of the loop, rather than selecting a few points along the loop, 
thus maximizing the amount of information taken from the coronal image. 

We address the question of validity of such a reconstruction, given the fact that the 
coronal field is probably non-linear force-free and that the superposition of linear force-free 
fields obtained for each individual loop would not, in general, be a force- free fiel d. To do so. 



we pe rform a series of tests on non-linear analytic force-free fields, described in iLow fc Lou 



(jl990l ). and as model loops we utihze projections of field lines on the photosphere. We 
compare several results of the method with the original field: 3D shapes of loops, local twist 
(coronal a), distribution of twist in the photosphere and the strength of magnetic field. 

We also find that the best-fit line may not be a global minimum of d{a,h), but one 
in a particular part of parameter space. We developed an algorithm that aids in locating 
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the appropriate region of the parameter space. We followed it manually, but it could be 
automated. The algorithm seems to improve the results on Low & Lou fields, but it probably 
does not describe every possible feature of the parameter space resulting from every single 
arbitrary magnetic field. The algorithm could be improved and expanded, based on further 
research involving other non-linear fields and real solar data. 

Based on Table [1] we draw the conclusions (i) that on the trial fields the twist is recon- 
structed with mean absolute deviation of at most 15% of the range of photospheric twist, 
(ii) that heights of the loops are reconstructed with mean absolute deviation of at most 
5% of the range of trial heights, and (iii) that the magnitude of non-potential contribution 
to photospheric field is reconstructed with mean absolute deviation of at most 10% of the 
maximal value. 

As shown in Fig. [TSl there is a typical understimation of twist when performing this 
procedure. Based on the experiments with Low &; Lou fields we conclude, that typically 

cereal ~ ^•'^'^OL found' 

We also demonstrate how this method can be apphed to real solar data, by doing the 
reconstruction based on the data from SOHO/MDI magnetograms and Hinode/XRT X-Ray 
images. The resulting field lines visually match the observed loops, have reasonable heights 
and self-consistent amoun t of twist of the magii itude that agrees with existing measurements 
of twist in active regions (jBurnette et al. 2004 ) . 
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9. Appendix, on Low & Lou fields. 



Low fc Loul (jl990l ) constructed a class of non-linear force-free magnetic fields beginning 



with an axi-symmetric field. An axi-symmetric, divergence-free magnetic field can be written, 
in all generality, as 



B = X V0 + QV(j) 



1 (IdA^ dA^ 
r sm u \r ou or 



(5) 



where A{r, 6) is the flux function and the azimuthal component is B^{r, 9) = Q{r, 6)/rsm 9. 
The poloidal components of the force- free condition, V x B = ctB, are satisfied only if Q 
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and a are each functions of the flux function alone 



Q{r,e) = Q\A{r,e)\ , a[Air,9)] = ^ . (6) 



The azimuthal component of the force-free condition 

r sm V ■ — 7)— = —aQ—— = — - , (7) 

is known as the Grad-Shafranov equation for flux function A. 

The Grad-Shafranov equation contains one free function Q'^{A) for which Low and Lou 
took a particular form. We generalize their choice to 

Q'(A) = a" , (8) 

where a and n are free constants. The absolute value signs, absent from the original formu- 
lation, are introduced here so that is real, and non-negative, even where the flux function 
is negative. Equation ([7]) can then be made homogeneous in r by proposing a solution 

A{r,e) = P(cos^)r-" , (9) 

for a still-unknown function P{fi)- Using this in the Grad-Shafranov equation, and defining 
fj, = cos 9, leads to the non-linear equation 

(1 - /i')^ + n{n + 1)P = -a^ (l + -) |Pp/" P , (10) 

for the unknown function P(/i). 

Eq. fllOj) has real solutions for any a and any n > —2; n need not be an integer. Boundary 



conditions, similar to those of ILow fc Loul (Il990l ). are -P'(l) = —10 and -P(l) = 0, so the 



solution is regular along the positive z axis. The solution will be regular along the negative 
z axis only when the solution satisfies the additional condition P(— 1) = 0. For a given value 
of n this condition will be satisfied only for certain choices of the eigenvalue a. (Integer 
choices of n always have one eigenvalue, a = 0, for which P(/i) is sin^ 9 times the Legendre 
polynomial of order n.) 

The final magnetic field, defined for 2; > 0, is constructed by rotating the axi-symmetric 
field by angle $ about the y axis and translating it downward a distance i. For rotation angles 
$ < 7r/2 and non- vanishing displacement, i > 0, the origin and what had been the negative 
z axis lie in z < 0, outside our domain. This means the regularity condition, P(— 1) = 0, is 
not needed to assure a regular magnetic field. We therefore make no restriction on P(— 1) 
and consider both n and a to be free parameters. 
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The function Q{A), required for the final field of eq. ([5]), is found from the square root 
of eq. ([8]). When the solution P{fi) changes sign there can be more than one choice of square 
root for which Q{A) is a continuous real functionjf] We focus on two such choices, we call 
signed and unsigned 

Qu{A) = aA\A\'/^ , (11) 
Qs{A) = a|A|i+i/" , (12) 

both of which satisfy eq. (jS]) when A{r, fi) = P{fi)/r"' is real. Taking the derivative of these 
functions gives the twist parameters for the two cases 

au{A) = + l^r^" , (13) 

as{A) = a (^1 + ^) |A|i/"sgn(A) . (14) 

The unsigned case has a single sense of twist determined by the sign of the constant a; the 
signed case has both senses provided A changes sign. 
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Loop # a, arcsec ^ hmax, arcsec 






0.012 


12.9 


1 


0.016 


0.0 


2 


0.014 


23.8 


3 


0.026 


9.5 


4 


0.008 


35.4 


5 


0.024 


33.9 


6 


0.020 


14.3 


7 


0.014 


12.7 


8 


0.010 


31.5 




0.010 


24.5 


10 


0.012 


27.4 


11 


0.010 


3.3 



Table 2: The results of the a-h fit to solar data. Here hmax is maximal height (not height at 
mid-point, as before in the text). Loops that don't seem to give a good fit to the data are in italic. 




Fig. 19. — Out of 12 loops that we selected, three, namely, loops 4, 5, 9 don't seem to give a good 
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2007-02-1 OTl 1 :22;02,j2SZ I loop § 2007-02- 1 0T1 1 ;22:02.J28Z I loop ^ 1 2007-02- 1 OTl 1 ;22;02.J28Z I loop § 2 




Fig. 20. — Parameter spaces for loops shown in Fig. [T9l 



